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Abstract. We provide a detailed theoretical and conceptual study of a planned 
experiment to excite Rydberg states of ions trapped in a Paul trap. The ultimate 
goal is to exploit the strong state dependent interactions between Rydberg ions to 
implement quantum information processing protocols and to simulate the dynamics of 
strongly interacting spin systems. We highlight the promises of this approach when 
combining the high degree of control and readout of quantum states in trapped ion 
crystals with the novel and fast gate schemes based on interacting giant Rydberg atomic 
dipole moments. We discuss anticipated theoretical and experimental challenges on the 
way towards its realization. 
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1. Introduction 



Coherent Rydberg physics with ensembles of ultra cold atoms is a flourishing and rapidly 
growing field that currently receives broad attention across area boundaries. The reason 
for this interest is rooted in the remarkable properties of highly excited atoms which 
offer a new approach for the study of strongly correlated many-body physics and the 
implementation of quantum information processing protocols [I] . Rydberg atoms, where 
the outer valence electron is excited to high-lying electronic states with large principal 
quantum number n, interact strongly via dipolar or van-der-Waals interactions, which 
can lead to interaction energies of several MHz over distances of several micrometers. 

The hallmark and key mechanism underlying coherent Rydberg physics is the so- 
called dipole blockade. This blockade hinders two atoms from being simultaneously 
excited when localized within a certain blockade radius: The interaction between two 
nearby atoms is here so strong that the doubly excited state is shifted out of the 
laser resonance |2]^4] . In dense gases this drastically influences the excitation behavior 
when atoms are photo-excited to Rydberg states. This effect has been experimentally 
demonstrated both in dense atomic ensembles [5j[6] as well as for two atoms stored 
in nearby optical dipole traps [7j[8], where the latter setup has enabled the first 
realization of entangling two-qubit gate operations for neutral atoms (9 10 . This 
remarkable experimental progress has stimulated a number of theoretical works: these 
include the exploration of novel quantum phases 11 12 , the study of the coherent 



real time dynamics of strongly driven Rydberg gases 13-15 , the development of new 



schemes for the efficient generation of entanglement 16 -20 , and the development of a 



quantum simulation architecture for open-system dynamics of complex many-body spin 
based on Rydberg atoms stored in optical lattices [22] or magnetic trap 



models 21 



arrays 



23 



In recent work 24 some of the present authors suggested the use of Rydberg 



excitations in ion traps to simulate many-body spin systems and to implement quantum 



information processing protocols. At present laser cooled and trapped ion crystals 25 



provide one of the most precisely controllable many body quantum systems available 

Ions can be addressed individually by laser light performing 

The 



in the laboratory 26 , 27 



28 



operations on the internal electronic structure representing a spin or qubit 
last decade has seen various experimental realizations of entangling quantum gates 
with a fidelity better than 99 percent 29 , the demonstration of simple quantum 



algorithms, entanglement swapping, teleportation 30 31 and recently the generation 



of entanglement with up to 14 ions 32 . Moreover, during the last years, the potential 
power of this ion trap quantum computer to emulate the physics of interacting many- 
body quantum spin-models [33] has been demonstrated in a series of remarkable recent 
experiments 



34 35 



Also, by combining coherent quantum gates and optical pumping 
as a dissipative element, the building blocks of an open-system simulator with trapped 
ions 



36] have been realized in the laboratory [371. 



Present experiments follow the paradigm of entangling qubits represented by 
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internal electronic or spin degrees of freedom of trapped ions in a linear chain via 



common vibrational modes of the trapped ions representing a phonon data bus 38 , 39 



For a linear (ID) ion trap the mode spectrum for large crystal sizes becomes increasingly 
complex, and scalability of ion trap quantum computers is achieved by shuttling ions 



around 40-42 and by going to 2D trap arrays 43,44. In contrast, gates based 



on Rydberg interactions do not rely on common vibrational modes, and thus the 
ionic Rydberg gate may provide an interesting alternative approach to scalability 
Additionally, Rydberg physics promises long-ranged and strong interactions, such that 
the interaction dynamics takes place on a nanosecond time scales, which is much faster 
than the external motion of the ion string. In Ref. 24 this idea was explored for 
the first time in some detail, and in particular also the dynamics of a single Rydberg 
ion in the presence of the electric trapping fields was described as a composite object 
consisting of a highly excited valence electron and a core formed by the nucleus and 
closed inner electronic shells. In the parameter regime of interest the electronic wave 
function extends over distances larger than the localization length of the ionic core, but 
is still considerably smaller than the inter-ionic distances of a few micrometers. Also it 
was proposed to enhance and tune the interaction strength of ions in Rydberg states by 



employing microwave dressing techniques 24 . In the context of quantum simulation 



of many-body spin models, it was argued that the resulting interionic interactions 
would allow one to realize spin chain dynamics with coherent excitation transfer on 
a nanosecond time scale. Likewise, these state-dependent long-range interactions might 
be explored for realizing fast two-qubit quantum gates, potentially several orders of 
magnitude faster than used entangling operations via the common vibrational modes. 
Thus the ability to excite and manipulate Rydberg excitations of trapped ions could 
provide a new interesting route for implementing quantum information processing 
protocols and the study of many-body physics in ion traps. 



1.1. Overview 

The purpose of this work is to undertake a quantitative case study of the feasibility to 
perform coherent excitation of Rydberg states in a crystal of 40 Ca + ions held in a Paul 
trap. Ultimately, this will lead toward the creation of the coherent state-dependent 
interaction of an ion crystal as sketched in Fig. [TJ We approach this task from three 
angles: 

(i) In Sec. |2]we discuss in detail the envisioned ion trap setup. For demonstrating the 
coherent excitation of single ions to Rydberg states we intend to create a trap with 
many ions configured in a linear string at the center of the trap or optionally an 
elongated three dimensional crystal. Here, the trapping of many ions will enhance 
the probability to excite ions and will simplify the task of aligning the narrow focus 
of the vacuum ultraviolet (VUV) laser to the ions. For the experimental study 
of Rydberg interactions at a later stage this trap will be operated by trapping a 
few ions in a string configuration with small distances between the individual ions. 
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Level structure 




Figure 1. Envisioned setup. Cold 40 Ca + ions confined in a Paul trap form a linear 
(a) or three-dimensional (b) crystal. Each ion may encode a qubit in which the logical 
states |0) and |1) are stored in the ground state 3Si/ 2 and the long lived 3-D 5 / 2 state, 
respectively. Two- or multi-qubit gate operations are realized by exciting a Rydberg 
state \r) of selected ions. Once excited these ions interact over long distances and 
on a fast time-scale, independently of the vibrational mode structure. The coherent 
excitation to the Rydberg state is achieved by a single VUV photon provided by tightly 
focused lasers. 



We show that it is possible to reach both settings with a standard linear Paul 
trap that provides free access for all laser beams and the observation of the ion 
fluorescence. Furthermore we explain the detection of single Rydberg excitations 
using an electron shelving scheme. 

Sec. [3] is dedicated to a detailed discussion of the laser setup that will be used 
to coherently excite ionic Rydberg states. In our envisioned scheme Rydberg nP- 
states are excited from the metastable 3D 5 / 2 state. One of the key challenges here 
is the required short wavelength close to 123 nm. This energy gap can in principle 



be overcome by using multi-photon excitation schemes 45 . We propose a different 



approach where the coherent excitation of Rydberg states is achieved by a single 
photon. We provide a study of how such a laser source is realized and how it is 
integrated into the ion trap setup. The key tool for such a continuous- wave (cw) 
VUV laser beam is four-wave mixing (FWM) in mercury vapor with fundamental 
laser powers of several 100 mW. This is already an established technique for 
production of a coherent Lyman-a beam tuned to the 1S-2P transition in hydrogen 
at 121.56 nm 46,47 . By choosing fundamental wavelengths close to resonances of 



mercury the FWM efficiency can be strongly enhanced. A laser source around 
123 nm with power of a few fiW should be feasible by using a triple resonant FWM 



scheme 48 . 
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(iii) In the theory part of this paper (Sec. |4j) we provide a thorough theoretical analysis 
of the quantum dynamics of a single trapped 40 Ca + ion that is excited to a 
Rydberg nP-manifold. We calculate the electronic energies and estimate the Rabi 
frequency for the laser excitation of Rydberg states. Our investigation reveals that 
the inhomogeneous field forming the Paul trap effectuates an inevitable coupling 
between the internal and external degrees of freedom of the ion. We calculate 
trapping field induced electronic level shifts as a function of the principal quantum 
number and quantify the expected change of the trapping potential for Rydberg 
states. Furthermore we analyze transitions among states within the nP-manifold 
that are induced by the oscillating trapping field. Our results show that the 
likelihood of such transitions strongly depends on the principal quantum number 
and on the phase of the radio- frequency field of the Paul trap. 



2. Ion trap 

In this section we discuss the setup of the ion trap which is designed for two categories of 
experiments with Rydberg ions. In the first step, we will investigate single ion excitation 
physics of 40 Ca + to Rydberg states in the time-dependent potential of the ion trap. The 
adaptation of spectroscopic methods for the Rydberg states, the measurement of the 
line width of the excitation resonance and the selection of the optimal Rydberg level 
for further experiments are specific goals here, as well as the detection of wave packet 
dynamics in the dynamic potential of the ion trap which is theoretically investigated in 
Sec. |4| Specifically suited for these experiments is a long string consisting of many ions 
as shown in Fig. [T^i, or even an elongated three-dimensional ion crystal loaded into the 
trap. This geometry in conjunction with a proper alignment of the laser will allow to 
increase the overall probability to excite Rydberg ions. 

The next step will be to study interaction effects between two ions when excited to 
Rydberg states. For these experiments we require ideally an inter-ion distance below the 
Rydberg blockade radius [3], which is accomplished by increasing the trap frequencies. 
In order to get similar Rabi frequencies for all ions, the ions should be configured in 
a linear string of a total length that is much shorter than the Rayleigh length of the 
focused VUV laser. The design of the trap allows to reach both settings with the same 
trap geometry. 



2.1. Experimental setting 

We will trap 40 Ca + ions in a linear Paul trap consisting of four cylindrical rods and 
two endcaps. The trap is shown together with the VUV laser setup in Fig. [2] The 
diameter of the rods is 2.5 mm while the distance r from the trap center to the rods is 
1.1 mm. This ratio between the diameter and ro was chosen to optimize the quadrupole 
approximation for the electric field of the trap 49 , thus enabling us to store three- 
dimensional elongated crystals. The distance between the endcaps of the trap is 10 mm. 
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four-wave mixing region separation interaction region 




Figure 2. VUV laser setup, beam guidance and excitation region. Lower 
part: Schematic of the fundamental laser systems. (SHG: second harmonic generation; 
LBO,BBO: nonlinear crystals; BSO: beam shaping optics; Ti:Sa: titamsapphire laser). 
Upper part: Vacuum system. The overlapped fundamental beams are focused into the 
mercury cell where four-wave mixing takes place. Separation of the VUV beam from 
the fundamentals is achieved using the dispersion of a MgF2 lens. A second MgF2 lens 
both focuses the VUV beam into the ion crystal and seals the low vacuum region (at 
~ 10~ 4 Pa) from the ultrahigh vacuum (UHV) chamber (at ~ 1CP 9 Pa) which contains 
the Paul trap. The VUV beam passes holes in the endcap electrodes of the trap and 
is monitored by a photomultiplier tube (PMT). 



Calculations 50 yield for an RF amplitude of U ac = 500 V at a frequency of £7 = 
27rxl5MHz, a stability parameter of q ~ 0.2 and a radial confinement u r of about 
27rxlMHz. To achieve an overlap of the beam with as many ions as possible, the 
VUV laser is aligned along the axis of the trap and propagates through holes in the 
endcaps. For the initial alignment of the VUV laser, we will aim for an elongated three- 
dimensional crystal with a few concentric shells around the axis. Following Ref. [51] , 
the density no of the ions in such a crystal can be estimated as (eoU ac ) / (Mfl 2 ^) where 
M is the mass of the ion. For the above trap parameters the distance between two ions 
is then approximately 11 «m. Comparing this value to the waist of the VUV beam at 
the focus of only 1.5 «m and a Rayleigh length of 50 /im one finds that only the central 
string of ions will contribute to the Rydberg yield once the beam is properly aligned, 
hence we will switch to a linear configuration of ions for the experiments with single 
Rydberg ions. 

The maximum number N of trapped ions in a single linear string configuration 



depends on the ratio between axial and radial trap frequencies 52 . The critical ratio 



is approximated by (u z /u r ) 2 = 2.94 x N , so the linear string configuration is stable 
if the radial frequency is much higher than the axial one. Simulations of the axial 




Figure 3. Dependence of the probability to excite a Rydberg ion on the trap 
settings. The left plot (a) shows the effective Rydberg excitation enhancement rate 
for a linear ion string compared to one single ion in the focus of the beam. Assuming 
the trap is loaded with the maximum number of ions for a linear chain configuration the 
excitation rate will increase for higher radial and lower axial trap frequencies. Within 
the given parameter region we can expect an enhancement by about a factor of 20. The 
right plot (b) shows the relation between the trap frequencies and the minimal distance 
between two ions. Distances down to 4 fim are possible with realistic parameters for 
the voltages and the RF frequency and the given trap geometry. 



trapping potential [53] result in uj z of about 27rx80kHz for a dc voltage of V& c = 90 V 
which leads, together with the radial trap frequency of 27rxlMHz, to a linear string 
consisting of 30 ions with a total length of 250 /zm (25). Summing over the relative 
excitation probabilities of the individual ions I z /Io, where I z is the laser intensity at the 
place of the ion and Iq the intensity in the focus, and assuming that the VUV beam is 
perfectly aligned to the ion string (see Fig. this setting allows to increase the total 
excitation probability to the Rydberg state by more than a factor of 15 as compared to 
spectroscopy with a single ion in the focus of the VUV beam. The detection scheme for 
the spectroscopy of Rydberg levels is outlined in Sect. |3.4| and we point out that the use 
of a large linear ion crystal is helpful since the overall excitation probability is increased 
and the alignment of the laser beam at 123 nm with its very tight focus becomes less 
critical. 

For later experiments investigating the dipole blockade between two Rydberg ions, 
the ions should be configured in a short linear chain on the trap axis. The distance 
between two ions has to be smaller than the blockade radius which is estimated to be 
of the order of 5 fxm [24]. The minimal distance between two ions depends on the axial 
confinement and the number of trapped ions and can be calculated [25] as 

/ e 2 \ 1/3 2.081 
min \^e Mul) X iV0-559- W 

If we increase the parameters to a radio frequency of 27rx30MHz with an amplitude 
of 2000 V and an endcap voltage of U^c = 1600 V, we get trapping frequencies of u r 
= 27rx2.2MHz and u z = 27rx400kHz. This would lead to an ion crystal with N = 
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12 ions and a minimal inter- ion distance of 4.1 /mi. The different settings for the ion 
trap are illustrated in Fig. [3] showing (a) the increase in total excitation probability to 
the Rydberg state for different axial and radial trap frequencies and (b) the minimal 
distance between ions for experiments exploring the Rydberg blockade. For a further 
increase of the trap frequencies in order to get to distances of 3 /mi and less, and in a 



later stage of the experiment, one might employ an advanced segmented micro trap 54 



3. Laser system 

The light required for the Rydberg excitation from the 3D 5 / 2 state to a nPj state 
lies in the VUV region between 131 nm (n = 10) and 122 nm (n — > oo). Four- wave 
sum frequency mixing (FWM) is an established technique to produce cw light in this 



frequency regime 47 [55 57]. A good candidate as nonlinear medium is mercury vapor 
which provides suitable transitions enhancing the four-wave mixing efficiency In this 
section we will first discuss the principle of threefold one-photon resonant four-wave 
mixing which strongly enhances the mixing efficiency. In the second part the setup of 
the laser source will be presented and in a third part we estimate the laser power for 
different ?>D^/ 2 -nPj transitions. 



7 1 S 



6 1 S 




Figure 4. Four- wave mixing in mercury vapor. The UV beam at 254 nm is close 
to the 6 1 5-6 3 P resonance, the blue beam at 408 nm establishes together with the UV 
beam the two-photon resonance (G 1 ^-? 1 ^). The green laser frequency determines the 
resulting wavelength of the four- wave mixed beam in the VUV. 



3.1. Principle of VUV generation 

Relevant levels of mercury and frequencies for the four-wave mixing process are shown 
in Fig. [4j The UV beam at 254 nm is red-detuned to the 6 1 S-6 3 P resonance. The 
blue beam at 408 nm establishes the two-photon resonance with the 6 1 S-7 l S transition 
together with the UV beam. The third beam in the green wavelength region determines 
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the output wavelength in the VUV due to four-wave mixing. The resulting power 



through four-wave sum frequency mixing is given by 58 

Here is the angular frequency and Pi the power of the zth beam, N the density of the 
nonlinear medium, xi the nonlinear susceptibility per atom and Ak a = {k^ — k\ — ki — 
&3V-V0 the wavevector mismatch per atom. G(bNoAk a ) the phasematching function for 
phase matching by the density of the nonlinear medium giving a maximum value at 
bN$Ak a = —4, where b is the beam 6-parameter given by two times of the Rayleigh 
length. For efficient four-wave mixing the only free parameter, besides increasing the 
fundamental intensities by using higher powers and tighter focusing, is the nonlinear 
susceptibility per atom x^ ■ This can be enhanced by using fundamental wavelengths 
close to resonances of mercury. Especially a two-photon resonance is essential for 



cw four-wave mixing 59 . Additionally we use two one-photon resonances to further 
increase the four-wave mixing efficiency. Therefore the UV wavelength is set close to 
the 6 1 S-6 3 P resonance and the VUV wavelength is set close to the 10 1 P and 11 X P 
resonance, respectively. This triple resonant scheme enhances the four-wave mixing up 
to 4 orders of magnitude and establishes a continuous VUV laser source in the /iW 
power range (48] . 



3.2. Laser Setup 

A schematic of the setup for the VUV production is shown in Fig. |2j The three 
fundamental beams are shaped by pairs of spherical and cylindrical lenses. The beams 
are overlapped at dichroic mirrors and focused into the mercury cell using a fused silica 
lens with a focal length of 15 cm. The mercury cell can be heated up to 240 °C providing 
a mercury vapor density of up to iV = 1.1 x 10 24 m~ 3 . Outside the focus region cooled 
baffles are used to avoid condensation of mercury on the optics. The produced VUV 
light is separated from the fundamental light by dispersion at a MgF 2 lens (/ = 13 cm 
at 123 nm and / = 20 cm at 254 nm) and a tiny mirror in the focus of the fundamental 
light. A second MgF2 lens (/ = 12.5 cm at 123 nm) separates the low vacuum region 
from the high vacuum region (trap) and focuses the VUV beam into the 40 Ca + trap 
with a Rayleigh length of 50 fim and a beam waist of 1.5 /im. The power in the VUV 
can be monitored by a photomultiplier behind the trap. 

A schematic of the fundamental laser system is shown in the lower part of Fig. [2] 
The beam at 254 nm is produced by a frequency-quadrupled Yb:YAG disc laser (ELS, 
VersaDisk 1030-50). Frequency-quadrupling is done with two resonant enhancement 
cavities, the first one using a lithium triborate crystal (LBO) as nonlinear medium, the 
second one using a /3-barium borate crystal (BBO). From 2W of infrared (IR) light at 
1015 nm we get up to 200 mW of UV radiation. This system is in principle capable 



to produce up to 750 mW of UV light, for details see 60 . The second fundamental 



beam at 408 nm is produced by a frequency-doubled titanium:sapphire laser (Coherent, 
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899-21), pumped by a frequency doubled Nd:YVC>4 laser (Coherent, V10). The external 
frequency-doubling cavity uses a LBO crystal. From 1.5 W of IR light at 816 nm we get 
up to 500 mW of blue light. 

For the desired VUV wavelength of 122 nm to 123 nm the third fundamental beam 
must be in the wavelength region 555 nm to 580 nm. Frequency doubling IR light 
produced by a Yb-fiber amplifier seems a feasible laser source for powers up to several 



watts in the green 61 62 . A frequency stabilized laser diode can be used as a seed 
laser. By splitting the amplification into two separated stages (a low power stage and 
a high power stage) we expect low amplified spontaneous emission and high damage 
threshold of the Yb-fiber. Frequency doubling in an external resonator with a noncritical 
phasematched LBO as nonlinear crystal should achieve a conversation efficiency of more 
than 60%. Frequency tuning of the green beam can be done by tuning the oscillator 
frequency and matching the LBO temperature. 

3.3. VUV Power Estimate and Linewidth 

For efficient four-wave mixing in the 123 nm region two resonances in mercury are 
relevant: 6 1 S-10 1 P at 122.04 nm and 6 1 S-1 l l P at 123.223 nm wavelength. Triplet 
states of mercury have a several orders of magnitude smaller dipole moment to the 
ground state and therefore a much smaller influence on four-wave mixing. Figure [5] 
shows a calculation of the four-wave mixing efficiency near the two singlet states. On- 
resonance the VUV radiation generated ist absorbed by the mercury vapor. In the 
proximity of a mercury resonance, however, the efficiency is much increased. The 10 1 P 
state in mercury yields an efficiency of about 50 /iW/W 3 for the transition to the 24P 3 / 2 
state. And even an efficiency of 150/iW/W 3 is achievable by using the 11 X P state of 
mercury which corresponds to the 67P3/2 state of 40 Ca + . For the calculation we took 
account of absorption of the VUV light in mercury vapor and imperfect fundamental 
beam qualities [63). With moderate fundamental powers of 200 mW (UV), 300 mW 
(blue) and 3 W (green) we would get between 9 /iW (n = 24) and 27 juW (n = 67). The 
corresponding Rabi frequencies are estimated in Sec. |4.4[ We expect saturation effects 
at several juW due to pumping of mercury out of the ground level but a /xW VUV source 
appears realistic with this setup. 

The laser linewidth of the VUV source depends on the emission bandwidths of the 
three fundamental laser systems. In the present state these laser systems are free- 



running. In an earlier experiment on Lyman-a excitation of atomic hydrogen 46 
linewidths very close to the natural linewidth of the atomic transition have been 
observed. This gives an upper limit of 10 MHz laser linewidth at 121,6 nm. For the 
planned experiment narrower linewidths are desired. We thus will frequency-stabilize 
the fundamental laser systems. A laser-linewidth in the VUV in the MHz range and 
below seems feasible. 
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Figure 5. The FWM efficiency near the mercury resonances. Calculated 
efficiency of FWM close to the ll 1 P (left) and 10 1 P (right) state of mercury. The dots 
mark positions of D 5 / 2 ~ n P3/2 transition wavelength. 



3.4- Ion laser cooling, excitation and detection of single Rydberg ions 

In order to excite the ions to Rydberg states they have to be prepared by Doppler- 
cooling and subsequently optical pumped into the 3D 5 / 2 state. The wavelengths of the 
required transitions for cooling and pumping are well suited for diode laser sources. 



Another advantage of the 40 Ca + level structure is the very long lifetime of 1.2s 64 of 
the intermediate 3_D 5 / 2 state which is the lower state in the Rydberg excitation process. 
Fig. [6] shows all relevant levels and transitions. For loading the trap with ions, we will 
employ photo ionization, with diode lasers near 423 nm and 375 nm or optionally by 
using a multi-photon ionization with light at 532 nm from a pulsed Nd:YAG laser. 

We propose a sequence of laser interactions, similar to standard ion trap quantum 
computing schemes, for the Rydberg excitation: First the ions are (i) cooled on the 
45*1/2 to 4Pi/ 2 transition near 397 nm close to the Doppler limit. Laser beams near 
866 nm and 854 nm provide for optical pumping out of the metastable D-levels. The 
emitted fluorescence near 397 nm is imaged on the chip of an EMCCD camera. In a 
next step (ii) we optically pump all the ions into one of the Zeeman states of the 5i/ 2 
ground state, switch off the laser at 854 nm and (iii) excite all ions in the crystal to 
the 3-D5/2 level, which exhibits a 1.2 s life time. For the excitation we will rely on the 



efficient rapid adiabatic passage pulse 65 . In an initial stage of the experiment one 



may incoherently optical pump into the desired metastable level via the excitation to 
the 4P 3 / 2 state at 393 nm followed by decay to the 3_D 5 / 2 level with means of a UV light 
emitting diode. Even when illuminated by resonant radiations near 397 nm and 866 nm, 
no fluorescence is emitted now. After the preparation of the ion crystal in the 3D^/2 level, 
all laser sources are switched off. The Rydberg state is then excited (iv) with a pulse 
of frequency v Rydherg near a wavelength of 123 nm. If an ion was excited successfully 



to the Rydberg state here, it will decay within a few 10 to 100 /xs 66 predominantly 
to the 45*1/2 ground state. Thus, after the excitation pulse, the laser beams at 397 nm 
and 866 nm may be switched on again (v) and the emitted fluorescence of the decayed 
ion can be detected with an EMCCD camera. Under typical operation conditions, we 
are able to discriminate between the on and the off state of fluorescence within 2 ms. 
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Figure 6. Excitation and detection scheme for 40 Ca ions. The ions are first 
pumped into the metastable D 5 / 2 state, from where they are excited to a Rydberg 
state by a pulse from the VUV laser. From the Rydberg state they will decay into the 
ground state. When the two lasers at 397 nm and 866 nm are switched on after the 
VUV pulse, the fluorescence of the ground state ion can be measured. 

The VUV laser near 123 nm pulse may optionally be triggered on the zero crossing 
of the radio frequency drive at frequency Q such that the ions are excited under the 
same electric field conditions in the trap. This may be important to study effects on 
the Rydberg electron wave packet in the time-dependent potential, theoretically studied 
in Sec. |4.5[ The entire scheme allows for a near to 100% discrimination of individual 
single Rydberg excitation events within a few ms and provides a spatial resolution of 
the diffraction-limited image of about 2 /im, much better than the inter-ion distances, 
such that we can determine the location and also the correlations of excitations in the 
ion crystal. At the end of the sequence a final optical pumping step (vi) with laser 
light at 854 nm, 866 nm and 397 nm brings ions in the whole crystal back to the S1/2 
ground state. This entire sequence will take about 5 to 10 ms and will be repeated as 
the laser frequency v Rydherg is tuned over the relevant frequency range of the Rydberg 
levels. Equally, one may vary the pulse length or apply Ramsey pulses in the Rydberg 
excitation step. 

4. Single ion theory 

This section is dedicated to a thorough theoretical investigation of the dynamics of a 
single trapped 40 Ca + ion. We will derive the underlying Hamiltonian, investigate the 
energy level structure of highly-excited Rydberg states and provide a detailed account of 
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the coupled electron-ion dynamics. In our study we will particularly focus on Rydberg 
nPj states which are laser-excited from the metastable 3D 5 / 2 state. We estimate 
the dipole-matrix element for this transition and calculate expected shifts of the trap 
frequency of the ion when excited to Rydberg states. Finally, we analyse possible state 
changing transitions which occur due to the inevitable coupling of internal and external 
dynamics of the ion. Simple scaling laws of the above-mentioned effects with respect to 
the principal quantum number n will be derived. 

4-1- Rydberg ion in a linear Paul trap 

Here we consider a 40 Ca + ion trapped in a linear Paul trap. The Hamiltonian governing 
the dynamics of the valence electron and the ionic core can be approximated by the 
two-body Hamiltonian 

H = H + + V (\ r * ~ R iD + V U T * - R i) + 2e$(Ri, t) - e$(r e , £)(3) 
2M 2m e 

where Pi (p e ) and M (m e ) are the momentum and mass of the ion (the electron) and 
e is the electronic charge. Ri (r c ) is the coordinate of the ion (the electron) and the 
potential V(|r c — Ri|) is an angular momentum dependent model potential which we 
will discuss in more detail in the subsequent subsection. We consider furthermore the 
spin-orbit coupling which is accounted for by the operator 

with ais ~ 1/137 being the fine structure constant. Here 1 and s are the orbital angular 
momentum and the spin operator of electron, respectively. The last two terms of 
Hamiltonian ^ emerge from the coupling of the charges to the electric potential of 
the linear Paul trap. This potential is composed of a static and a radio frequency (RF) 
component reading 

$(r, t) = a cos ttt(x 2 - y 2 ) - (3(x 2 + y 2 - 2z 2 ). (5) 

Here a and (5 are the electric field gradients of the RF and the static field, respectively, 
and Q is the corresponding RF. The combination of these two fields permits the 
confinement of positively charged particles in the trap center. 

It is convenient to treat the system in the center-of-mass frame, in which the 
Hamiltonian ^ becomes 

H = H l + H e + H le . (6) 

Here Hi, H e and H e are the free Hamiltonian of the ion, the Hamiltonian of the electron 
including electron-trap coupling and the coupling between electron and ion, respectively. 
All of them are explicitly written as 

Hi =^ + e<l>(R,t), (7) 

H e =^+V(\T\)+VUr)-e^(T,t), (8) 
H le = -2e[a cos Vtt(Xx - Yy) - f3(Xx + Yy- 2Zz)\ . (9) 
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where corrections due to the finite nuclear mass are neglected. The electron-ion coupling 
Hi e results in an additional potential for ion. Our calculation shows that for energetically 
low-lying states, roughly n < 10 (see Sec. 4.6) this additional potential is far smaller 
than the effective trapping potential and thus can be neglected. The static field 
in conjunction with the rapidly oscillating RF field then gives rise to an effectively 



static (ponderomotive) harmonic potential for the 40 Ca + ion 67 and the corresponding 
Hamiltonian of the external motion can be approximated by 

P 2 M 

ZJW z P=X,Y,Z 



Here oo r = u x = oo Y = \j2[(ea/MQ) 2 - e(3/M] and u z = 2^e(3/M are the 
transverse and longitudinal trap frequencies, respectively. Using the trap parameters 
a = l.OxlO 9 V/m 2 ,/3 = 1.0 x 10 7 V/m 2 and Q = 2ttx25 MHz, one obtains u r = 2ttx3.27 
MHz and Uz — 2-7T x 1.56 MHz. These values will be used throughout the theory 
part. They are slightly higher than the ones discussed in Sec. [2] and will make the 
influence of the trap on the electronic dynamics more prominent. As mentioned before, 
neglecting the electronic structure when considering the trapping is only valid for low- 
lying electronic states in which the 40 Ca + ion can be regarded as a singly charged 
point-like particle. In the Rydberg states in which we are interested here, the spatial 
extension of the electronic wavefunction can easily become larger than the oscillator 
length of the trapped ion. To accurately describe ionic Rydberg states one thus has to 
account for the coupled electronic and external dynamics. 

Our strategy for doing this is to first solve the electronic Hamiltonian excluding 
the electric trapping potential. This is achieved by diagonalizing the electronic model 



potential 68 which yields accurate eigenenergies and corresponding wave functions. 
These states will then constitute a basis in which we can expand the ion-electron 
coupling. The problem of finding the eigenstates and eigenenergies as well as the time- 
evolution of the system is then reduced to the solution of coupled matrix equations. 



4-2. Calculation of Rydberg energies and wave functions 

We will now briefly outline the calculation of the field-free electronic wave functions. 
In this section we use atomic units throughout unless mentioned otherwise. The single 
electron in the outer shell of the 40 Ca + ion moves in the Coulomb potential of the 
nucleus that is screened by the inner electron shells. In general this is a complicated 
many-electron system but since the inner electrons are forming a closed shell it is 
possible to approximately reduce it to a two-body problem. Here the valence electron is 
orbiting in an angular-momentum dependent (i.e. dependent on the quantum number I) 
model potential that approximates the interaction with the nucleus and the inner shell 
electrons. Such potential can be parameterized as [68] 

Vi(r) = [2 + (Z - 2) exp {-a lA r) + a i>2 r exp(-a /j3 r)] 

-^{l-exp[-(r/r0 6 ]} (11) 
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where Z is the nuclear charge, a.\j and r; are parameters depending on angular 
momentum, and a cp is the experimental dipole polarizability of the doubly charged 
ion. All these quantities are tabulated in Ref. [68). The spin-orbit interaction is then 
described by the Hamiltonian 



v£ o) (r) 



c&l-sdV^r) 



2r 



dr 



a 



1 -2 



V,(r) 



(12) 



where the last factor is introduced to regularize the wave function at the coordinate 
origin. The model potential conserves the total angular momentum j = 1 + s, which 



permits an expansion of the electronic wave function into the spinor basis 69 



(13) 



where 4>i,j(r) is the radial wave function and \l,j,rrij) denotes the two-component spinor 
basis vector 



\ l 3 m j) = Yl ('» V 2 5 m ^ m sli> m j)l^ m ;)l s ^)- 

m s =±l/2 



(14) 



Here {I, 1/2, mi, rrij — mi\j,mj) are the Clebsch-Gordan coefficients and rrij (mi) is the 
quantum number of the projection of the angular momentum j (I) onto the quantization 



axis. Substituting Eq. ( 13 ) into the corresponding Hamiltonian and integrating over the 



angular variables, we obtain the radial Schrodinger equation 

hi,j(j>i,j(r) = e ttj (j)i tj (r), (j)ij(0) = 0, 
with the radial Hamiltonian 



hi. 



1(1 + 1) 



2r 2 



+ V,(r) 



^T\r) 



(15) 



(16) 



We solve this equation by discretizing space and diagonalizing the corresponding matrix. 
The electronic energies obtained from this calculation are in good agreement with the 
experimental ones 70 . Hence we can be confident that the obtained electronic wave 



functions form a solid foundation for the upcoming analysis of the Rydberg states of 
the trapped ion. 



4-3. Bare electronic level structure 

Before treating the coupled electron-ion problem we will first discuss the energies of 
the bare electronic states obtained from solving the eigenvalue problem of the radial 



Hamiltonian (16). We will focus on Rydberg P-states as they are most relevant to the 



envisioned experimental excitation procedure. The Rydberg P-states of 40 Ca + are split 
into two manifolds with total angular momentum j = 1/2 and j = 3/2. Due to the 
quantum defect these states are furthermore energetically well isolated from adjacent 



levels. According to Ref. 71 one can parameterize the Rydberg energy, including spin- 



orbit coupling, according to 



E, 



nlj 



2(n - 5 lj f 



(ai s Z c ) 2 ( 1 3 
(n-tfy) + 1/2 4(n-fy). 



(17) 
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Figure 7. Bare electronic energies of the nPj Rydberg states. The energy of 
the states P\ii and P 3 / 2 in the field- free case is measured with respect to the state 
3D* function of the principal quantum number n. The inset shows the fine- 

structure splitting between nPi/2 and nP 3 / 2 - Note the different energy units of the 
inset. 



Here Z c is the net charge of the ionic residue seen by the valence electron. The second 
term in the brackets gives rise to the fine structure splitting due to the spin-orbit 
interaction. The quantum defects of the nPj are found by fitting the data with Eq. 
(17), #1,1/2 = 1-4396 and #1,3/2 = 1.4358. We have calculated them with respect to the 



state 3D5/2, whose energy can be found in Ref. [72] and in the NIST database. The 
huge energy gap, typically a few PHz (1 PHz corresponds to a laser wave length of 
Xl ~ 299.8 nm), requires a VUV laser for exciting the transition 3D 5 / 2 — > nPj. The 
energy difference (fine-structure splitting) between nPi/ 2 and nP 3 / 2 , which is shown in 
the inset in Fig. [TJ is well-approximated by AE n ij w j^f^i 71 ~~ ^j)~ 3 which neglects 
differences in the quantum defect for different principal quantum number. The fine- 
structure splitting varies from a few tens to hundreds of GHz when n < 20, to a few 
GHz when 20 < n < 50 and reaches sub GHz when n > 50. This has to be taken into 
account when nPj-states shall be excited selectively from the states 3D 5 / 2 a s the laser 
linewidth has to be considerably smaller than the energy gap between nP\/ 2 and nP 3 / 2 . 
It is instructive to compare these energy scales with other energies of the system. 



We use the parameters of Sec. 4.1 i.e., the RF field has the frequency f2 = 2n x 25 MHz 
and the oscillation frequencies of the ion in the effective harmonic trap are ujz ~ 27r x 1.56 
MHz and u r ~ 2ir x 3.27 MHz with the parameters provided above. Thus, neither the 
RF field nor the vibrational motion can couple the P1/2 to the P3/2 manifold resonantly 
if n is not too large, i.e. AE n ij 3> Q 3> u r ~ uz- 
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4-4- Laser excitation of Rydberg states 

The Rabi frequency for the envisioned laser excitation of Rydberg states from the state 
3D 5 / 2 is determined by the dipole matrix element (3D 5 / 2 | r | nP^^. In order to obtain 
an estimate for this matrix element we calculate the radial wave function of the 3D5/2 



with the model potential (11). We find for principal quantum numbers > 18 that 

;3D 5/2 (3/2) \z\nP 3/2 (3/2)) « 0.383 x a ^ L588 (1* 



with a being Bohr's radius. To estimate the expected Rabi frequency we use the laser 
power estimates given in Sec. |3.3 A beam waist of 1.5 /im and a laser power of 9 /iW 



(n = 24) and 27/iW (n = 67) lead to a Rabi frequency of 2ir x 5.2 MHz (n = 24) and 
2tt x 1.8 MHz (n = 67), respectively, when the ion is placed in the laser focus. We 
have to point out that the present calculation 3D 5 / 2 wave function does not fully take 
into account correlations of the inner electrons and the valence electron. Hence, the 
presented value can only be seen as a first estimate and a more rigorous calculation will 
be needed for more accurate results. On the one hand the numerically obtained energy 
of 3D5/2 agrees well (within 0.04%) with the experimental value given in Ref. 70 . On 



the other hand, we find significant deviations (up to 8 %) when comparing numerically 



calculated oscillator strengths with experimental values tabulated in Ref. 73 . 

This concludes our discussion of the level structure of 40 Ca + in the field- free case. 
Before eventually turning to the discussion of the coupled electronic and ionic dynamics, 
we will investigate in the following two subsections the electron-trap and the electron-ion 
interaction in order to gain a further understanding of the system. 



4-5. Electron-trap interaction 

To calculate the interaction of the valence electron with the electric fields of the Paul 
trap one has in principle to expand the entire Hamiltonian ^ in the electronic basis 
and then subsequently diagonalize it or propagate the desired initial wave function. 
This is particularly difficult when dealing with Rydberg states due to a high density 
of states which will lead to a high dimensional matrix. To circumvent this problem 
we exploit - as earlier mentioned - that the nP-manifold is well isolated from other 
states. This permits the derivation of an effective Hamiltonian that acts only in 
this low-dimensional electronic subspace but takes into account the couplings to other 
electronic states perturbatively. The derivation of this Hamiltonian is based on the 
van Vleck transformation [74| . This procedure greatly reduces the dimension of the 
matrices involved and therefore simplifies the calculation and analysis of the problem 
considerably. The effective electronic Hamiltonian becomes six-dimensional. 

We now turn to the calculation of the electronic level shifts due to the electron-trap 
interaction in Eq. ([8]). The interaction Hamiltonian can be written explicitly in the 
spherical coordinates and reads 

H et = -A^e(3r%%9, ( p)-2J^ea r 2 co S nt[Yi(e^) + Yi 2 (9, ( l ) )}.(19) 



Rydberg excitation of trapped cold ions: A detailed case study 



18 



This means we are actually calculating the electronic level-shift for a situation in which 
the ion is located perfectly at the trap center. In this subsection we consider the time 
t merely as a parameter. This quasi-static approximation seems reasonable since the 
frequency of the RF field is typically much smaller than the frequency associated with 
the motion of the Rydberg electron. At a later stage we will however take a full account 
of the time-dependence when describing the coupled electron-ion dynamics. 

As indicated in the introductory paragraph of this section the level shifts caused by 
H et are much smaller than the distance between the field-free electronic eigenenergies 
E n ij. We can thus calculate them by using second order perturbation theory. 
The corresponding van Vleck transformation yields the following effective interaction 
Hamiltonian in the reduced nP,-subspace: 

Here the sum over v = {n, I, j, m,j} excludes all states in the nP,-subspace of interest and 
jj(X) mv = (iJ) m \H et \il) v }. From the selection rules imposed by the structure of Eq. (19) 
[the orbital angular momentum of two coupled states has to differ by ±2] it is evident 



that the matrix elements H^ mu will not couple any electronic state within the nP, 



3 



subspace directly. The coupling however can happen indirectly through intermediate 
states as shown by Eq. (20). Here our numerical calculations show that the diagonal 
matrix elements of Eq. (20) are always larger than the off-diagonal elements at any 
moment. For example, at a = 10 9 V/m 2 and (3 = 10 7 V/m 2 , the ratio of the diagonal 
elements over the off-diagonal elements of the state 35P varies in range of 4 to 200. 
Therefore the electronic level shifts are mostly determined by these diagonal entries. 

The numerical results for the diagonal matrix elements, as functions of the electric 
field gradients a and (3 for the manifold 35P, are depicted in Fig. 18} When VLt = ir/2, 
the terms involving a actually vanish. Here the energy shifts H ( ^ nm are determined 
solely by (3, i.e. the static electric field of the trap. For the range of principal quantum 
numbers shown, a variation of (3 causes a shift of the electronic energies which is of 
the order of a few hundred kHz. However, since usually a ^> (3 the RF field is causing 
much stronger level shifts which for the data shown can be up to a few MHz when 
| cos(fit)| = 1. This time-dependence of the electronic energies has to be accounted for 
when the Rydberg states of a trapped ion are laser excited from low-lying states. 

For all times one finds that the states of the 35P are shifted to a lower energy 
compared to the bare electronic state, arise from level repulsion from states above nPj. 
Hence the electric field mostly causes mixing of the nPj with the nearby manifold of 
P-states. The energy shifts H^ mm depend rather sensitively on the principal quantum 
number n. From Eq. (19) and Eq. ([20 ) one finds that they scale according to n 11 , due to 
(r 2 ) ~ n A and AE n ij ~ n~ 3 . In Fig. 8 we show numerical results (for | cosf2t| = 1) that 
indeed confirm the scaling. Energy shifts up to several hundred MHz, can be found when 
n > 60. Treating these high lying states accurately, however, requires to go beyond the 
second order in the perturbation expansion. Moreover, such high principal quantum 
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Figure 8. Electron-trap interaction under the assumption of a quasi-static 
RF field. Here we plot the diagonal electronic energy shift u^ mm as a function 
of the principal quantum number n for fit = kn with integer k (log-log scale), (a) - 

(b) Dependence of the electron-trap interaction on a and /3 at different time t. (a) 
fit = kir; here the RF field reaches its maximal amplitude and causes the strongest 
level shifts, (b). fit = kn + tt/2; here terms involving a vanish and only the static 
field component of the Paul trap (with magnitude (3) determines the energy shifts. 

(c) . Magnification of the energy shifts in the region n = 30 to n — 35. In (c) linear 
scale is used. The parameters for which the data is calculated are fl — 2tt x 25 MHz, 
a = 10 9 V/m 2 and /3 = 10 7 V/m 2 . 



numbers necessitate the consideration of field ionization a discussion of which can be 



found in Ref. 24 



4-6. Electron-ion coupling 

We will now focus on the actual electron-ion interaction which is effectuated by the 
inhomogeneous nature of the electric fields forming the Paul trap. Similar couplings are 
also observed in context of Rydberg atoms that are trapped by means of inhomogeneous 



magnetic fields 75 , 76 . For trapped ions the corresponding coupling Hamiltonian is H\ e 
[Eq. (§]. 

To treat it we will again make use of the van Vleck transformation which yields the 
following matrix elements of the effective coupling Hamiltonian in the nPj-subspace: 



(2)mm' 
Ie 



RJU R 



where R is the position of the ion and 

v -' mm' / j 1 *"e *"e 



+ 



(21) 



(22) 



2 (-Em — E v ) 2(E m i — E v ) 

Here we have used the abbreviations H^ mu = (ip m \H e \ip u ) and R c = 2e[(— a cos fit + 
f3)x, (a cos fit + {3)y, —2(3z}. When the electron is in the low-lying states the energy 
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shifts caused by Eq. (21) are negligible in comparison to the typical spacing between 



vibrational energies of the ion. This is due to the fact that the coupling is roughly 
proportional to the size of the electronic orbit since it depends on the relative coordinate 
r of the ion-electron system. The situation is changed when the electron is excited to 
Rydberg states. Here the large orbits of the electron, (r) 



a n 2 , leads to substantial 



change of the magnitude of Eq. (21 ) 




18 36 
Principal quantum number n 



Figure 9. Estimated trap frequency Awj of the additional potential induced 
by the electron-ion coupling in Rydberg states. The dependence of Acox on 
the principal quantum number n is shown (log-log scale). The curves shown scale 
according to ~ n 7/2 . Here Q = 2tt x 25 MHz, a = 10 9 V/m 2 and /3 = 10 7 V/m 2 . 



The coupling caused by Eq. (21) can affect the motion of the ion in two possible 



ways: first, the effective time averaged potential that governs the external ionic motion 
is not separable anymore, because the off-diagonal matrix elements of U mm / couple the 
ionic vibrations in different directions. Second, the diagonal elements of U mm ' will 
lead to an additional trapping potential for the ion that is superimposed to the already 
present harmonic confinement. 

To quantify the strength of the additional trapping potential that is effectuated 
by the electron-ion coupling, we calculate the time-averaged energy shift AEx" 1 of the 
electronic energy as a function of the displacement in X direction from the coordinate 



center. Using Eq. (21) the result is 

2 (2a 2 + 4/3 



2\ y{2)mmj^2 



(23) 



where y^ mm [ s a parameter that depends only on the electronic states. It scales as 
The harmonic frequency of this additional potential can be estimated in accordance with 



~ n 



7 
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Eq. (ph to yield 



2e W + 4/3 2 )Vg )mm 

N m • (24) 



This result does not take into account modifications of the trap frequency that arise from 
cross-couplings to the oscillatory motion in the Y and Z direction and is thus merely 



an estimate. Note that Aw x ~ a/« 2 + 2(3 2 ~ a and Aux ~ yVi^ mm ~ n 7 / 2 . This 
scaling indicates that the oscillation frequencies of a trapped ion can be substantially 
modified when the electron is excited to Rydberg states. In Fig. [9] we show some 
numerical results obtained for 40 Ca + . One observes that for very high lying Rydberg 
states the trap frequency Aux of the additional potential becomes even comparable 
to Ux itself. This indicates that a description of the external ionic motion through 
harmonic oscillators might break down entirely. 



4-7. Adiabatic potential surfaces 

In the previous two sections, the individual electron-trap and electron-ion interactions 
were investigated. This analysis provided us with basic insights into the coupled 
ion-electron dynamics. We are now in position to study the complete interaction 
Hamiltonian in Eq. ^ and Eq. (|9]). Such treatment poses difficulties in particular 
when dealing with the oscillating RF field. The conventional approximations usually 
made in case of such oscillating fields, e.g. the rotating wave approximation, are hard 
to carry out in the present problem, as quite a few different timescales are involved (see 
section 4.3). 




Figure 10. Adiabatic potential surface of the 35P3/2. The potential surface is 
plotted along the X — £lf-plane. The dashed curves are lines along which degeneracies 
occur. Here n = 2tt x 25 MHz, a = 10 9 V/m 2 and f3 = 10 7 V/m 2 . 
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To make progress we employ once again the approximation of a quasi-static electric 



field which was already used in section 4.5, i.e., we calculate the eigenvalues of the sum 



of the Hamiltonians ^ and ^ while taking time as a parameter. As result we obtain a 
four-dimensional (three spatial and one temporal axis) adiabatic potential surfaces. In 



Fig. 10 we show a cut (along X and fit) through the surfaces that belong to the states 
35P3/2- We do not show the contribution of the ionic center of mass Hamiltonian (|7j) 
which just gives rise to a state-independent potential. 

We observe two pairs of degenerate potential surfaces. This degeneracy is typical 
for Hamiltonians representing half-inter spin particles and a time-reversal invariant 
Hamiltonian. It is called Kramer's degeneracy [77|. In addition all four potential surfaces 
become completely degenerate (in the X — fit-plane) when cos fit = /3/a. Note that the 
complete degeneracy can only occur in the X — Y plane at particular moments in time, 
namely when the static and RF fields cancel. These degeneracies induce non-adiabatic 
transitions between neighbouring electronic states, which forbid the application of a 
single-surface Born-Oppenheimer approximation for the description of the coupled ion- 
electron system. One thus has to solve the complete quantum problem consisting of the 
three external ionic degrees of freedom plus the six-dimensional electronic state space 
of the nPj-manifold. 

4-8. State- changing transitions in the nPj-manifold 



The degeneracies shown in Fig. 10 can in principle lead to non-adiabatic transitions 
among the electronic states of the nPj-manifold. Such state-changing transitions are 
certainly not desirable in particular when ionic Rydberg states are excited for the 
purpose of implementing quantum information processing protocols. We will now study 
the severity of this effect in more detail by analyzing the temporal evolution of a wave 
packet in which the ion is initially prepared in the state 

|*(0)> = |000)|35P 3/2 (3/2)), (25) 

where 1 000) is the ground state of the effective harmonic oscillator Hamiltonian ( JToj ) (to 
characterize the vibrational state we use the notation \nxnynz) where n& corresponds 
to the number of vibrational quanta in the fc-th direction). The state |^(0)) can be 
approximately created by a fast 7r-pulse on the transition 3D 5 / 2 (5/2) — > 35P 3 / 2 (3/2) 
applied to an ion that is cooled to the external ground state. We will study the evolution 
under the action of the Hamiltonian ^ where we replace Hi by the effective Hamiltonian 



(10). This approach neglects the micromotion. We are interested in the population of 
states other than |^(0)) function of time. 

In Fig. 11 we show the corresponding plot calculated for the parameters fl = 
2n x 25MHz, a = 10 9 V/m 2 and /3 = 10 7 V/m 2 . The time interval shown corresponds 
to twice the period of the transverse external ionic motion 4n/u r . Apart from the 
numerically exact result we also present the loss of population calculated from time- 
dependent perturbation theory. Both curves are in excellent agreement, showing only 
a marginal "loss" of population from the initial state which oscillates periodically at a 
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Figure 11. Population dynamics. (a) - (b) The initial electronic state is 
|35P 3/2 (3/2)). (a) Population loss from the initial state |000)|35P 3/2 (3/2)). (b) In the 
course of the time evolution population is transferred to the states |200) |35P 3 / 2 (3/2)) 
(red) and |020)|35P 3 / 2 (3/2)) (black). The electronic state is virtually unchanged, (c) 
- (d) The initial electronic state is |50P 3 / 2 (3/2)). (c) Population loss from the initial 
state |000)|50P 3/2 (3/2)). (d) The time evolution of dominant modes |200) |50P 3/2 (3/2)) 
(red) and 1 020) |50P 3 / 2 (3/2)) (black). During the time interval shown, no significant 
transition to the states nPjjtj . In all the panels solid/dashed curves correspond to 
numerical simulation/time-dependent perturbation theory. The parameters used in 
the simulation are, Q = 2n x 25MHz, a = 10 9 V/m 2 and (3 = 10 7 V/m 2 . 



frequency of about 2tc x 6.7 MHz. The maximum loss, which is about 0.01%, is reached 
at about 70 ns. The fast oscillations visible on top of the curve are due to the RF 
field. The absence of loss for this particular set of parameters is partially because all 
the transition matrix elements are small. For instance, the maximal off-diagonal matrix 
elements of the electron-ion interaction is here about 2tt x 400 kHz, which is several 



times smaller than the diagonal matrix elements in Eq. (21). Moreover, the RF field is 



far off resonant with all the transition energies, as shown in section 4.3 



A further analysis of the wave packet propagation shows that essentially only the 

are populated in the course of the 
(b). This demonstrates 



11 



states |200) 35P 3/2 (3/2)) and |020) 35P 3/2 (3/2) 
time evolution. The corresponding data is presented in Fig. 
that virtually no transitions among electronic states are taking place and that solely the 
population of the external trap levels is slightly redistributed. 

Though population transfer is small for the state |35P 3 / 2 ), it can become significant 
for higher Rydberg states. The matrix elements for transitions between vibrational 
modes are proportional to U mm i and scale proportional to n 7 . As a result, the maximal 
population of the excited states can be estimated to be proportional to n 1A as long 
as perturbation theory applies. As an example, we calculated the loss dynamics for 
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the state |50P3/2(3/2)). The corresponding data is shown in Fig. 11 (c) - (d). The 
fully numerical calculations show that the loss of the population from the initial state 
is about 143 times larger than the one of the state |35P3/2(3/2)), which is in a good 
agreement with the scaling. Moreover at higher Rydberg states the dramatic increase 
of the transition matrix elements and the effective trap frequency give rise to higher 
oscillation frequencies, which is clearly shown in the dynamics of the dominant modes, 



also depicted in Fig. 11 (c) - (d). Consequently, the redistribution of the population of 



the external modes becomes more important at the higher Rydberg states. 
5. Conclusion and Outlook 

So far our study revealed no conceptual obstacles that speak against the feasibility to 
coherently excite Rydberg states of trapped 40 Ca + ions. The calculation of the spectral 
data presented in this work heavily relies on the accuracy of the model potential that 
was used for calculating the Rydberg wave functions [68]. First experiments will thus 
allow us to assess the accuracy of these data. These results will be fed back into 
the calculations for further improvement and refinement of the theoretical predictions. 
Most importantly, experiments will measure accurately the transition frequencies and 
the dipole matrix element for the transition 3P5/2 — > nP 3 / 2 for which we could only 



provide an estimate in Sec. 4.4 



This, together with accurate spectroscopic data will direct the next step in which we 
intend to perform a detailed study of the interaction properties between ionic Rydberg 
states. As in the single ion case the consequences of the coupling between the internal 
and the external ion dynamics caused by the inhomogeneous trapping field will have to 
be understood. We expect that, unlike in the case of neutral atoms, the van-der-Waals 
interaction among excited ions will be too weak to achieve large (~MHz) interaction 
energies. To overcome this problem we intend to study the possibility to "dress" Rydberg 



with microwave photons. In Ref. 24 it was shown that this is a promising route towards 
strongly interacting spin-models and fast quantum gates. Quantum state tomography 
and quantum process tomography have been well established when using the low lying 
qubits states S1/2 and D 5 / 2 in 40 Ca + . This technique will allow to characterize the 
entanglement operations by the Rydberg excitation. In the longer run, one might 
address single ions with the tightly focused VUV beams in the setting indicated by 
Fig. [T] (b), analogue to the experiments with neutral atoms trapped in focused laser 
beams (7|||. 

Beyond that, the envisioned experimental setup might be potentially used to study 
novel hybrid quantum states in ion traps, where one highly excited electron is shared 



among two or more core ions 78 . This would pave the way toward an entirely new 



regime in trapped ion physics. Finally, Rydberg states of ions could be coupled to 



microwave resonators and strip line cavities, as discussed for atoms in Ref. 79 , but now 



as hybrid solid state - trapped ion systems in the spirit of Ref. 80 . 
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